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Abstract 

Dark matter particles gravitationally trapped inside the Sun may annihilate into Standard Model 
particles, producing a flux of neutrinos. The prospects of detecting these neutrinos in future 
multi-kt neutrino detectors designed for other physics searches are explored here. We study the 
capabilities of a 34/100 kt liquid argon detector and a 100 kt magnetized iron calorimeter detector. 
These detectors are expected to determine the energy and the direction of the incoming neutrino 
with unprecedented precision allowing for tests of the dark matter nature at very low dark matter 
masses, in the range of 10 — 25 GeV. By suppressing the atmospheric background with angular 
cuts, these techniques would be sensitive to dark matter - nucleon spin-dependent cross sections at 
the fb level, reaching down to a few ab for the most favorable annihilation channels and detector 
technology. 
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I. INTRODUCTION 



A plethora of cosmological and astrophysical measurements over the last decades has 
tested the validity of the Standard Model of Big Bang cosmology (a spatially flat Friedman- 
Robert son- Walker model) at an unprecedented level of precision. Current observations point 
to a flat universe in which the mass-energy content includes only 5 % of ordinary matter 
(baryons) while 22 % is a mysterious non-baryonic dark matter (DM) component, whose 
properties are largely unknown. If it can be identified with some relic, this should be a 
massive, neutral and sufficiently long-lived particle. A popular hypothesis, realized in many 
extensions of the Standard Model of particle physics (SM), including supersymmetry 1| 
or extra dimensional models as 

, is that the DM is dominated by a single species of 
weakly interacting massive particle (WIMP). Such WIMPs could have masses ranging from 
several GeV to dozens of TeV, see Ref. {4] for a review. Direct DM searches look for the recoil 
energy of target nuclei due to interactions with the DM particles. WIMP scattering on nuclei 
may happen via spin-dependent interactions (proportional to J(J + 1), rather than to the 
number of nucleons) or via spin-independent scattering (which grow with A 2 and therefore 
increase dramatically with the mass of the target nuclei). The spin-independent cross section 
dominates over spin-dependent scattering in current direct DM detection experiments, which 
use heavy atoms as targets, setting therefore very strong constraints on spin-independent 
interactions. For DM capture in the Sun, however, the distinction is not so strong as in 
direct detection searches since a large fraction of the Sun's mass is in the form of hydrogen, 
for which there is no enhancement of the spin-independent cross section. In the following, 
we will present our results in terms of spin-dependent DM scattering cross section. However, 
in practice, similar neutrino fluxes would be obtained with spin-independent cross section 
around two orders of magnitude smaller due to its enhancement for the heavier nuclei present 
in the Sun (see, e.g., Ref. |5|). Thus, all the results presented here for spin-dependent 
scattering cross sections can be roughly translated for the spin-independent case with an 
improvement of around two orders of magnitude. 

Although there is no confirmed observation of the passage of WIMPs through a detector, 
a long standing positive signal in the form of an annual modulation has been reported 
by the DM search experiment DAMA/LIBRA 6j. This annual modulation is expected in 
the scattering of DM particles off nuclei in the detector due to the Earth's annual motion 
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around the Sun, that modifies the relative velocity of the Earth and the DM halo. The 
DAMA/LIBRA collaboration has reported the observation of 13 such cycles, compatible 
with the expected phase due to the Earth's orbit, at 8.9cr. This observation is, however, 
not confirmed by other direct search experiments, such as Xenon and CDMS fl. 

Tension between the two datasets remains even when interpreting the DAMA/LIBRA results 
as scatterings of DM particles off the Na nuclei 1 , which would point to light DM masses 
~ 10 GeV with a spin-independent scattering cross section a ~ 0.2 fb close to the threshold 
of the Xenon and CDMS sensitivities. On the other hand, the CoGeNT experiment has 



also reported an observed excess of events over the expected background [10J. As in the 
DAMA/LIBRA case, the shape of the excess would also be consistent with scatterings of light 
DM ~ 10 GeV, although with somewhat smaller spin-independent scattering cross section 
a < 0.1 fb. The similarity of the DAMA/LIBRA and CoGeNT results is very intriguing, 
although there is some tension between the preferred value of the cross section and there is 
controversy on whether the two DM hints can be reconciled in the simple spin-independent 
elastic scattering picture, with analyses pointing to a combined allowed region around ~ 
7 GeV and a ~ 0.2 fb U| or disfavouring such an explanation [jjzj]. Even more tantalizing is 



an annual modulation, similar to 



13|. Furthermore, the CRESST-II 



the recent announcement of the CoGeNT collaboration o: 
the DAMA/LIBRA one, in their observed excess at 2.8cr 
experiment also observes an excess of 32 events to be compared with an expect background 
of 8.4 ± 1.4 in the O band fl^ . corresponding to the lowest energies. When interpreted as 
a DM hint, these events could be accommodated with a mass and cross section in the same 
region as that favoured by DAMA/LIBRA and CoGeNT 12j. While the accumulation of 
hints favouring light DM with masses of ~ 10 GeV and a spin-independent cross section 
a ~ 0.1 fb is tantalizing, this region is disfavoured by the null results of the Xenon and CDMS 
searches. However, as this region is close to the energy threshold, these results are sometimes 
questioned since small changes in the energy calibration or efficiencies could indeed reopen 
the window required for the DM interpretation of DAMA/LIBRA, CoGeNT and CRESST. 
On the other hand, the Xenon and CDMS experiments have a cleaner signal region, given 
the stronger background suppression achieved. Therefore, the DAMA and CoGeNT excesses 



1 The DAMA detector exploits very pure sodium iodide crystals and DM scatterings off the I nuclei would 
instead point to heavier DM masses. 
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could be due to a misinterpreted background component. A very similar situation occurs 
when interpreting the DAMA/LIBRA signal through a spin -dep endent cross section in the 



pb range and the null results of the PICASSO 



15|], COUPP [16| and SIMPLE [171 searches. 



Clearly, more experiments sensitive to these low DM mass regions are needed in order 
to settle the present tensions and shed light on the DM puzzle. In this context, DM indi- 
rect detection methods offer a very attractive complementary probe. They aim to detect 
the gamma-ray, positron, anti-proton and neutrino fluxes produced in DM annihilations in 
regions with dense DM accumulations, such as the Galactic center. In the case of neutrinos, 
the higher densities of DM expected to be gravitationally trapped in celestial bodies, such 

offer a very appealing alternative since, contrary 



as the Sun 



18 



19| or the Earth 



20 



to positrons or gamma-rays, neutrinos can escape the solar interior providing a snapshot of 
the interaction taking place. These indirect searches, provide a very complementary tool 
to direct detection probes since, unlike them, they will not only provide a measurement of 
the DM mass and scattering cross section, but also information on the different annihilation 
branching ratios, depending on the particular channel searched for. Bounds on the WIMP 
mass and cross section from the water Cerenkov (WC) neutrino detector Super-Kamiokande 
(Super-K) have been obtained for different annihilation channels, see Ref. 5|, |22|, |23| . A lot 



of effort has also been devoted to using neutrino telescopes to detect neutrinos from DM 



annihilation processes 



24j-|29[. However, these neutrino telescopes are mostly sensitive to 



and events, and their energy threshold is ~ 5 — 10 GeV in the muon energy (which 
is translated into ~ 15 — 20 GeV neutrino energy). Other neutrino detectors might be nec- 
essary to test the DM mass region indicated by DAMA/LIBRA and CoGENT via indirect 
searches. Even if Super-K is sensitive to low DM masses, the energy resolution achieved by 
the former detector is not as good as for some future facilities. Future neutrino detectors, 
such as WC, liquid argon (LAr), liquid scintillator, or magnetized iron neutrino detector 
(MIND) facilities, might be constructed for other physics measurements (such as neutrino 
oscillation studies or core collapse supernova neutrino physics) and could be exploited for 
indirect DM detection. 

As first proposed in Ref. 30j, future MIND-like detectors currently under development, 



such as India-based Neutrino Observatory (INO), could provide an excellent energy and 
angular resolution to neutrinos from WIMP annihilations in the Sun. In such detectors, 
the expected degeneracies among the DM mass, annihilation branching ratios, and elastic 
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cross section with nuclei could be resolved, shedding much light on the DM nature. Other 
detection techniques originally proposed and constructed for long baseline and supernova 
neutrino studies have also been examined in this context, such as future LAr detectors 
as well as current and future liquid scintillator detectors 32 ]. 
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Here we extend the work of Ref. [30J by considering both lower DM masses (~ 10 GeV) 
and additional annihilation channels to the bb and t~t + previously studied. We also consider 
not only future MIND detectors, but also proposed future LAr facilities, such as GLACIER. 
Both techniques are expected to achieve unprecedented energy and angular reconstruction 
capabilities at ~ GeV neutrino energies. 

The structure of this paper is as follows: Section [TT] describes the MIND and LAr detector 
technologies exploited here. The neutrino fluxes and event rates from WIMP annihilations 
in the Sun, as well as backgrounds to our signal from atmospheric neutrino events in these 
two detector facilities are presented in Sec. IIHI In Sec. IIVt we present the DM nature 
reconstruction exercise performed here by means of a Markov Chain Monte Carlo (MCMC) 
analysis for several light DM particle masses. We summarize and draw our conclusions in 
Sec. El 



II. LARGE NEUTRINO DETECTORS 



Large underground neutrino detectors 33M39I] have been very successful in observing neu 



trinos from several natural and artificial sources and have provided one of the few evidences 



we have of physics beyond the SM, as well fundamental results in astrophysics 
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Active research is underway worldwide to plan and build the ne xt-g ener at ion of large un 



derground neutrino detectors in the mass range of 10 2 to 10 3 kt 42j-|48|, which will play a 
crucial role in the upcoming era of precision neutrino science. The physics potential of such 
large detectors is extremely broad. One of the major aims is the search for proton decay 49], 



which provides a direct test of Grand Unified Theories. These large neutrino detectors are 



also essential for the future accelerator based long base 
to explore mixing and CP violation in the lepton sector 



ine neutrino oscillation experiments 
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45 2j. Applications in astrophysics 



include precise measurements of supernova neutrinos and the search for new sources of as- 



trophysical neutrinos 



43 



53[. In this context, the possibility of observin g th e neutrinos 



originated from the annihilation of DM particles in the center of the Sun 
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FIG. I: The angular resolutions for LArTPC and MIND considered in thi s study as a function of 
the true neutrino energy. These curves have been derived from 
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541, 



and to reconstruct the annihilation process so as to probe the nature of DM is particularly 
tantalizing and will be the main goal of our study. These detectors are expected to measure 
the energy and the direction of the initial neutrino with very high precision which allows to 
search for DM particles in the mass range of 5 — 50 GeV, which cannot be probed effectively 



with large neutrino telescopes such as IceCube 5 

Currently, four different technologies for large underground neutrino detection are being 
investigated worldwide: WC, unsegmented liquid scintillator, LAr time projection chamber 
(LArTPC), and MIND. Below we discuss in detail the LArTPC and MIND alternatives, 
which can provide the best performance for the searches under study. Indeed, the energy 
and angular resolution capabilities of WC and liquid scintillator detectors degrade rapidly 
beyond the quasielastic regime and thus makes them unsuited for the multi-GeV energies 
under consideration. The assumed characteristics of the detectors used in our numeric 
simulations are summarized in Tab. [H 
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Detector characteristics 


MIND 
(Only M ±) 


LArTPC 
(Both /z± & e ± ) 


Fiducial mass 


100 kt 


34/100 kt 


v energy threshold 


2 GeV 


2 GeV 


Detection efficiency (e) 


Fig. 6.21 (Top Panels) 
of Ref. [54] 


100 % for /i ± 
80 % for e ± 


v energy resolution (SE) 


0.15J3 


0.15^/1 GeV GeV 


v angular resolution (59) 


See Fig. U 


See Fig. U 


Bin size 


1 GeV 


1 GeV 



TABLE I: Detector characteristics used in our simulations. 



A. Liquid Argon Time Projection Chamber 

A LArTPC [56 1 is a modern version of the bubble chamber technology. The working 
principle of LArTPC is that by applying an electric field in highly pure argon, free electrons 
created by the passage of an ionizing particle can be drifted over large distances, C(m), 
without any distortion. This enables the possibility to obtain the three dimensional image 
of the particle track by just reading out the surface of the volume. The projections of 
the track can be used to reconstruct the three dimensional path of the particle with an 
accuracy of a few mm. The viability of the LArTPC technology has been tested by the 
extensive R&D program, developed by the ICARUS collaboration. Currently, the largest 



LArTPC ever built, is the ICARUS T600 module 



57 



59| with a mass of 600 1. Based on 



this technique, a very massive next-generation 100 kt detector (GLACIER) |43|, |48|, |60| has 



been proposed under LAGUNA 61] . The LBNE physics program in United States is also 



considering a 34 kt LArTPC as a possible detector candidate for DUSEL |51l . |62| . In this 
work, we study the prospects of a 34/100 kt LArTPC detector to search for DM WIMPs 
with 10 full years of data taking. 

Since this technology is still in its R&D phase, many uncertainties regarding its perfor- 
mance at very large mass scale exist. Thus, the detector characteristics that we consider 



here, are essentially based on the performance of the ICARUS T600 detector 



the ongoing estimates of the LAGUNA |6lj and LBNE [62 



57H59j and on 



631 ] collaborations. 
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In a LArTPC, the muon-electron misidentification probability is almost zero and it offers 
very high detection efficiency for muons and electrons in the GeV energy range, suppressing 
the neutral current 7r° backgrounds. In our study, we consider u e , u e , and charged 
current (CC) events with a neutrino energy threshold of 2 GeV. In principle, the technol- 
ogy can be extended down to 500 MeV, but these low energies are highly populated with 
atmospheric backgrounds. The highest DM mass that we consider in this study is 25 GeV 
and in this energy range one can expect that the muons and electrons produced via CC 
processes will be fully contained in these massive detectors. We assume an 80 % detection 
efficiency for electrons and positrons and for /i ± we consider a full 100 %. Electron tracks 
would remain fully contained up to very high energies enabling to probe DM masses up 
to 100 GeV. Since magnetization of this detector would be challenging, the charge of the 
particle is not measured and we combine the e~ and e + as well as the \x~ and fi + samples. 
We assume an energy resolution of SE(E) = 0.15^/1 GeV GeV for CC /i ± and e ± events. 
The atmospheric neutrino background can be suppressed to a very low level by selecting 
a narrow sky window of the size of the angular resolution of the detector centered around 
the Sun. The reconstruction of the incoming neutrino direction in a LAr detector has been 
estimated using the information coming from all the final state particles and the detector 
angular resolution as a function of incident neutrino energy has been shown in Fig. [H This 



information has been obtained from Fig. 4 of 



31] . Notice that by using information on both 



the leptonic and hadronic final states of the interaction, angular resolutions smaller than 
the typical scattering angle between the neutrino and the charged lepton can be achieved, 
something impossible if only the outgoing charged lepton is detected. We use the same 
angular resolution to estimate the atmospheric neutrino backgrounds for v e jv e and Vp.lv p. 
CC signal events. 



B. Magnetized Iron Calorimeter detector 

Iron calorimeters are made up of iron (steel) modules interspersed with sensitive elements 
in which charged particles deposit energy. Electrons and positrons cannot pass through iron 
and therefore these detectors are blind to v e /v e . This type of detectors is well suited to 
observe the long muon tracks and hence are capable of observing only Vp and Vp events. 
A magnetic field can be added to these detectors to distinguish Vp events from Vp. A 
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Magnetized Iron Neutrino Detector (MIND) 64|, |65j of large mass (50 - 100 kt) has been 
proposed for the precision study of neutrino oscillations of neutrinos from muon decays at 



a future neutrino factory facility 



52 



66 



for the INO facility 



70 



69(. A similar type of detector is being considered 



7jJ which is planned to have a total mass of 50 kt at startup, which 



might be later upgraded to 100 kt. The main aim of this detector is to measure atmospheric 
neutrinos 46[ and it may serve as a neutrino factory far detector at later stage 721 ] . 



In our study, we consider a 100 kt MIND detector with an energy dependent efficiency 
as given in Ref. 54J. For z/ M , it remains almost flat (~ 63%) in the neutrino energy range 
of 10 to 25 GeV while for z/^, it deteriorates from 76% to 71% going from 10 to 25 GeV. 
This detector can provide accurate measurements of the muon energy 54[. Moreover, due 
to the calorimetric nature of the detector, the energy of the hadronic shower can also be 
reconstructed j^, Q] which helps to reconstruct the energy of the interacting neutrino. We 
use an energy resolution of 5E(E) = 0.15E for CC fi^ events. The angular resolution of 
this detector is also quite good and mainly governed by the hadronic shower 55[. We use 
the information from both muons and hadronic showers 5J, |55[ to reconstruct the incoming 
neutrino direction. The neutrino angular resolution of this detector as a function of true 
neutrino energy has been depicted in Fig. [I] The angular resolution for is slightly better 
than v^. We use this information on angular resolution to suppress the atmospheric 
and backgrounds in a similar way to the one discussed above for the LArTPC, i.e. by 
discarding all events except those reconstructed within the solid angle subtended by the 
neutrino angular resolution centered around the Sun. Notice that the efficiency as well as 
the energy and angular resolutions of the MIND detector will depend on the number of 
active layers the track crosses and, therefore, on the direction of the incoming neutrino. 
While for a neutrino beam experiment the detector can be oriented along this direction this 
would be more challenging for neutrinos from the Sun. We have neglected this dependence 
in this first study since it will critically depend on the unknown orientation and latitude of 
the detector. 



III. NEUTRINO FLUXES AND EVENT RATES 

If DM particles can pair annihilate we can envision a wide variety of possible channels 
into lighter SM components. The rates for these annihilations will be greatly enhanced in 
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regions with larger density of DM particles. If the cross section between DM and nuclei 
is high enough, the collisions of DM with solar matter will allow them to lose energy and 
become gravitationally bound, accumulating in the solar gravitational potential. Thus, we 
can expect a larger number of DM annihilations in the Sun than in the surrounding space. 
Most of the annihilation products will be lost in scatterings and decays in the dense solar 
interior. Fortunately, neutrinos once again provide us with a unique tool to explore these 
secretive processes that could take place in the interior of the Sun. Indeed, neutrinos from 
DM annihilations could escape the solar medium unaffected and provide a picture of the 
processes taking place in the star's interior. If these neutrinos have energies in excess of GeV 
they will be clearly distinguishable from the MeV neutrinos produced in solar reactions. 

Many of the possible DM annihilation channels could indeed lead to such neutrino sig- 
nals 75]. Assuming that DM annihilations are flavour conserving and predominantly to two 



body final states, the bb, cc, t~t + or u a u a channels would lead to the required high energy 
neutrino probes of solar DM annihilations. For high enough DM masses, additional channels 
leading to high energy neutrino fluxes would also become kinematically allowed such as tt, 
W + W~, ZZ or even Higgs production. The energies of these neutrinos would make them 
unsuitable for the detector technologies under study and more adequate for large neutrino 
telescopes, since the associated charged leptons will not be fully contained in a small detector 
volume and the energy and angular resolution would be greatly degraded. We will therefore 
limit our study to DM masses up to 25 GeV. Other possible channels such as e~e + or 77 
would not lead to neutrino signals but could provide complementary probes in gamma ray 
or positron searches. 

The expected neutrino fluxes per annihilation at the detector from the different DM 



annihilation channels have been obtained using the WimpSim 



76| software. The fluxes for 



the "standard oscillations" with 13 = 0°, 9 23 = 45°, 9 12 = 33.2°, Am^ = 2.2 • 10~ 3 eV 2 
and Amy = 8.1 • 10~ 5 eV 2 were used. These fluxes must later be normalized with the 
annihilation rate in the Sun. In equilibrium this annihilation rate is one half of the capture 
rate, which is proportional to the DM-nucleon cross section and can be computed given 
the solar density profile, for which we used the standard BP2000 Solar Model (77)], a local 
DM density of Pdm = 0.3 GeV cm -3 , and a Maxwellian velocity distribution shifted by 
the solar velocity v® = 220 km s _1 . For the computation of the capture rate we follow l|, 
. Notice that large systematic astrophysical uncertainties affect the computation of the 
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FIG. II: Expected muon events in bins of 1 GeV width for the MIND (left) and GLACIER (right) 
100 kt detectors and for neutrinos from the annihilation of DM particles of 10 GeV (up) and 
25 GeV (down) mass. Ten full years of data taking for a 1 fb DM-nucleon spin-dependent cross 
section have been assumed. Since GLACIER has no charge identification both [i + and fjT events 
where added, while only \i~ events are displayed for MIND. 



capture rate and thus the overall normalization of the expected signal. We will not include 
this systematic error in the numerical results presented in Sec. UV] since their effect in the 
sensitivity to the annihilation cross section to the different channels is just a trivial rescaling 
by the corresponding amount. These fluxes were later convoluted with the cross section 
and detector efficiencies and energy resolutions summarized in Tab. [Jin order to obtain the 
expected number of event rates depicted in Fig. [Til The event rates have been computed for 10 
full years of data taking with a DM flux computed assuming a 1 fb cross section between DM 
and nucleons. For each channel a 100 % branching ratio is shown. For the MIND detector 
the muon neutrino events arising from the different annihilation channels are depicted. Since 



11 



MIND is a magnetized detector, muon antineutrinos can also be distinguished. The event 
rates for antineutrinos are very similar in shape but slightly lower than the neutrino ones, so 
we do not reproduce them again in a separate figure but will include them in the numerical 
analysis presented in Sec. IIVI On the other hand, GLACIER is not magnetized and neutrino 
and antineutrino events need to be summed. However, GLACIER is also suited to detect 
electron neutrino events. In Fig. [Til we show the muon neutrino and antineutrino events for 
the GLACIER detector. Electron events have a similar shape and we do not show them but 
they are also included in the analysis of Sec. IIVI As expected, the annihilation channels that 
lead to larger neutrino fluxes at the detector are the neutrino channels, with the hardest 
spectrum; the r channel with an intermediate spectrum; and the heavy quark channels (b 
and c) with a softer spectrum and somewhat smaller flux. The gluon channel (g), while also 
leading to neutrinos, provides an event rate too low to be observed and will be disregarded. 

Notice that all the neutrino channels provide a very similar spectrum (since the kinematics 
of the process is identical) and, moreover, neutrino oscillations ensure that almost equal 
amounts of all neutrino flavours will arrive at the detector. The and v T channels are 
indeed indistinguishable (since #23 was set to maximal mixing), and the v e annihilation 
channel only shows some deviation at low energies where the event rate is too low to be 
observed. Thus, when trying to reconstruct the DM annihilation channels in Sec. IIVI we 
will not distinguish between the different flavours but only consider a generic "neutrino 
channel" for which the case was chosen. Similarly, the annihilations into b and c quarks 
provide very similar spectra. As will be shown with a particular example in Sec. \1V\ the 
detector energy resolution is not enough to discriminate between these two channels and an 
effective "quark channel" q will be considered instead, for which the b case was chosen as 
representative. 

A common feature of all the channels is that the corresponding event rates will have an 
abrupt cutoff at the energy corresponding to the DM mass, providing a clean and unam- 
biguous measurement of this parameter. This will be true as long as either the r or neutrino 
signal is measured, given their harder spectrum. Thus, in the analysis of Sec. IIVI the DM 
mass will be fixed, since it is not degenerate with any other parameter. The free parameters 
that we will try to reconstruct will be the product of the DM-nucleon cross section times 
the branching ratio of DM annihilations to neutrinos, r and heavy quarks. 
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A. Atmospheric Neutrino Background 



The main source of background for the search of neutrinos from DM annihilations in the 
Sun will be the atmospheric neutrino events at the detector. Indeed, unlike solar neutrinos, 
the atmospheric neutrino energy range overlaps with that of the signal. This is an unavoid- 
able and potentially large background which can fortunately be reduced by exploiting the 
extremely good angular and energy resolution of the detectors under consideration. Indeed, 
while the signal neutrinos come from the direction of the Sun, the atmospheric neutrinos 
reach the detector from all directions and can be significantly reduced by accepting only 
the background events for which the reconstructed incoming neutrino direction corresponds 
with the solar position at the time of the event. 

In order to simulate the atmospheric neutrino background we take as representative fluxes 
those from Ref. 79[ for the Frejus site with mountain over detector at a solar maximum. 



These fluxes are then oscillated with the same oscillation parameters as the DM signal and 
convoluted with the cross section, detector efficiency and energy resolutions summarized in 
Tab. [H In order to estimate the background suppression achieved by the angular cut around 
the solar position, we first estimate the number of atmospheric events for a given energy 
integrating over the whole sky and then multiply the number of events by the solid angle 
subtended by a cone of the detector angular resolution at the corresponding energy divided 
by 4tt. Of course this is only an approximation and the atmospheric flux should only be 
integrated over the position of the Sun at a particular time. However, since this will depend 
on what is the actual geographical site selected for the detector, we regard this average 
over all the sky as a good approximation of the magnitude of the atmospheric neutrino 
background. We have considered only CC events in estimating the atmospheric neutrino 
background and in the numerical analysis we add a 10 % systematic uncertainty (fully 
correlated in all the energy bins) to this background. The contribution from neutral current 
(NC) events is negligible for both electron and muon samples due to the high background 
rejection capability (> 99.5 %) of the detectors under consideration. The atmospheric u T /u T 
fluxes produced via neutrino oscillations can give rise to tau leptons through CC interactions 
inside the detector. These taus can decay further into electrons and muons with a branching 
ratio of ~ 17 %. We have checked that these backgrounds are negligible. 

The black dashed lines in Fig. |TT] depict the estimated atmospheric neutrino background 
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for each detector after the angular cut. As can be seen, the background decreases fast with 
the neutrino energy and mainly affects the observability of the quark annihilation channels. 
Indeed, we have checked that increasing the background mainly affects the sensitivity to 
this channel, while the neutrino and r ones remain essentially unaffected. 



IV. RECONSTRUCTING THE DARK MATTER ANNIHILATION BRANCHING 
RATIOS 



As motivated in the previous section, we will scan and analyze the parameter space 
consisting of BR x cr {x = r, u, q) in order to determine the capabilities of the two detectors 
under study to discriminate between the different DM annihilation channels. This would 
provide invaluable information to shed light on the DM nature beyond the mass and cross 
section measurement that is available at direct searches and complementary to the other 
annihilation channels explored in indirect searches with gamma-ray satellites and telescopes. 
We will also explore the sensitivities that both detectors have to these parameters. The 



parameter space has been sampled using a MCMC based on the MonteCUBES [80| software. 
We perform this scan for simulated true values of both a = and BR x <r = 0.3 fb. The 
running time of the experiments were set to 10 years each with the efficiencies and energy 
resolutions summarized in Table [H For the MIND detector we separate the signal in muon 
neutrinos and antineutrinos, exploiting the magnetization of the detector to separate the two 
samples. For the GLACIER detector magnetization would be challenging and we therefore 
combine the neutrino and antineutrino samples. On the other hand, we exploit the capability 
of GLACIER of discerning the neutrino flavour and consider electron and muon neutrinos 
separately. In all cases, the threshold was set to 2 GeV and the events were distributed in 
bins of 1 GeV width. 

Since we are simultaneously scanning over all the parameters, we can easily see whether 
or not there are any correlations among the different branching ratios and what are the 
actual sensitivities to the parameters when marginalizing over the others. Regarding the 
DM mass, which could a priori present degeneracies with the other parameters under study, 
its determination is expected to be free of degeneracies with other parameters since it cor- 
responds to the cutoff of the neutrino spectra and will be easily detected for annihilations 
into neutrinos or taus. Thus, we will keep this parameter fixed in the present study for 



14 




BR b a [fb] 

FIG. Ill: The most probable regions at 68 %, 90 %, and 95 % posterior probability and a DM mass 
of 10 GeV projected onto the BR^o" — BR c <r plane. This simulation was made for GLACIER with 
a DM mass assuming a true BR^cr = 0.3 fb as the only source of neutrinos in DM annihilations. 
We can clearly see the linear degeneracy between the two parameters, motivating the study of only 
an effective quark channel BR^cr. 

simplicity. As we will show, in general, the sensitivities to the different branching ratios 
turn out to be fairly independent, since the sum will be bounded from above by the total 
event rate. 

For every scan in the MCMC, we generate four chains with a total number of 2 • 10 5 
samples per chain, check the variance within each chain and compare it with the variance 
of the overall sample in order to gauge the convergence of the chains. In all cases presented, 
we obtained a good convergence (R — 1 < 10 -3 ) 2 . The MCMC priors were assumed to be 
flat in the scanned parameters. 



A. Results 

In Fig. IHIl we illustrate the degeneracy between the DM annihilations into b and c quarks 
with a simulation for GLACIER with a DM mass of 10 GeV and BR^cr = 0.3 fb and keeping 
BRb<j and BR c a as free parameters in the fit. As we will see later, this corresponds to the 
most favourable setting and, from this figure, it is apparent that the two channels cannot 



2 See Ref. 8l| for the definition of this measure. 
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be distinguished, and thus it makes sense only to fit the combination BR g = BR;, + &BR C 
in our analysis, where k ~ 0.4 accounts for the smaller number of neutrinos resulting from 
each annihilation into cc. In order to simulate this channel we take the events expected from 
DM annihilation into bb. 

In Tab. |TTJ we present the results for our MCMC simulations for the case when there 
is no signal. The numbers correspond to the 90 % limit on the posterior probability and 
can be interpreted as the sensitivity reach of the experiment to each parameter. From the 
table, we can clearly see that the GLACIER detector would be significantly more sensitive 
to DM annihilations at low DM mass than the MIND one, having typical sensitivity of 
a few 10~ 2 fb for both BR r er and BR^a, where the corresponding MIND sensitivies are 
0(0.1 — 1) fb. There are several reasons why the MIND sensitivities are worse than the ones 
for GLACIER. First of all, the signal statistics are higher at GLACIER, since it considers 
both electrons and muons, and because of a higher assumed efficiency. In particular, the 
MIND efficiency is very bad at low energies, leading to the significantly worse bounds for the 
case of 10 GeV DM mass. Furthermore, the angular and energy resolutions of the LArTPC 
are better, leading to a better background rejection in these experiments. 

These sensitivity limits should be compared to the current bounds on the spin-dependent 
DM cross section. The most stringent direct bounds at the simulated masses are a < 400 fb 
at 10 GeV (PICASSO |lj) and a < 90 fb at 25 GeV (COUPP 16|). For the 10 GeV 
case, the present exclusion limit of PICASSO is also close to the lower part of the region 
required to explain the DAM A/LIBRA signal when interpreted as spin-dependent scatter- 
ings. Thus, depending on the dominating annihilation channel, the current bounds could be 
improved by up to five orders of magnitude in the most optimistic case of GLACIER, a DM 
mass of 10 GeV, and annihilations directly into neutrinos. However, it should be pointed 
out that DM is not required to have any branching ratio at all into channels that produce 
a neutrino flux, and thus direct detection experiments will always be complementary to 
indirect searches. However, if an annihilation branching ratio of at least ~ 10~ 3 to heavy 
quarks, ~ 10~ 4 to t + t~ or ~ 10~ 5 to neutrinos is present, the DM spin-dependent scattering 
interpretation of the DAM A/LIBRA signal can be probed at the GLACIER detector. Sim- 
ilarly, for a spin-independent scattering interpretation of DAMA/LIBRA, branching ratios 
of ~ 10 _1 to quarks, ~ 10~ 2 to t + t~ and ~ 10~ 3 at GLACIER could be probed. 

In comparison to the results obtained for Super-K [5], we find that the MIND experiment 
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would be of comparable sensitivity, while any of the LArTPC experiments studied would 
out-perform Super-K in most cases. However, it has to be borne in mind that the capabil- 
ities of a liquid argon detector have been significantly less studied and simulated than for 
magnetized iron calorimeters such as MIND, thus, the better performance of GLACIER is 
conditioned to achieving the optimistic assumptions listed in Tab. [I] Typically, the sensi- 
tivity to BRqCr is significantly worse due to the softer spectrum of this channel that makes 
it more contaminated by the atmospheric neutrino background. Indeed, we have verified 
that increasing the atmospheric background with less restrictive cuts greatly deteriorates 
the reach in BR^a while having a small impact in the measurement of the other two param- 
eters. Moreover, as can be seen from Fig. [Til a smaller flux of neutrinos per annihilation 
as compared to the other annihilation channels is obtained for the quarks. Conversely, the 
present bounds on this channel are also weaker and a MIND detector would actually pro- 
vide some improvement on this channel as compared to the Super-K bounds given its better 
angular resolution and, hence, background rejection. It is also interesting to note that the 
sensitivity of the GLACIER detector decreases with the DM mass assumed, as expected 
from the reduced neutrino flux (c.f. Fig. [TT|). On the other hand, the sensitivity of the MIND 
detector increases with the DM mass demonstrating the better optimization of this detector 
to high neutrino energies. 

In Fig. [IV] we show the preferred regions of the parameter space in the case of simulated 
BR-eO" = 0.3 fb. As expected from the sensitivities, the MIND detector is only capable of 
distinguishing BR„<7 from zero, while the signal for BR r cr is close to being detectable. On 
the other hand, GLACIER would be able to clearly identify the nature of the DM branching 
ratios into both rs and us, and being close to the sensitivity needed for detecting BK q a at this 
level. The reproduction of the parameters by GLACIER is remarkable and essentially free of 
correlations between the different branching ratios. There is a very small correlation between 
BR T and BR g , since both of these spectra are both softer than the BR„ one (c.f. Fig. ILTI) . 
but this correlation is almost negligible. Since the different branching ratios of DM could 
be singled out in an experiment like this, we would gain valuable information on the DM 
couplings to different particles, which in turn would be of large interest to constrain different 
models of DM. This would not be true in, e.g., direct detection experiments, which would 
only be sensitive to the total cross section on the target. It is important to remark that the 
discerning capability of the detectors between the different DM annihilation channels stems 
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Experiment DM mass 


BR T <7 [fb] BR u a [fb] BR g cr [fb] 


MIND (100 kt) 10 GeV 
25 GeV 


0.70 0.35 3.4 
0.34 0.15 1.7 


LArTPC (34 kt) 10 GeV 
25 GeV 


0.15 0.11 0.73 
0.16 0.10 0.21 


GLACIER (100 kt) 10 GeV 
25 GeV 


1.5 -10~ 2 6.4- 10" 3 0.25 
1.0 • 10~ 2 5.2 • 10~ 3 0.19 


Super-K data [5j 10 GeV 
25 GeV 


0.65 0.12 10 
0.45 0.19 5.0 



TABLE II: Sensitivities to the branching ratios multiplied by the capture cross section at 90 % 
posterior probability for different experiments and DM mass after 10 years data taking. The 
Super-K sensitivities are taken from Fig. 4 of Ref . . 



almost exclusively from their energy resolution and the ability to reconstruct the different 
spectral shapes shown in Fig. [TT1 On the other hand, the capability of the MIND detector 
to separate the neutrino and antineutrino samples does not help with this measurement 
since the corresponding antineutrino event shapes are extremely similar. Similarly, the 
information on the neutrino flavour GLACIER provides is largely irrelevant for the analysis 
and its ability to detect v e in addition to only provides an increase in statistics but no 
complementary information to solve degeneracies. The slight degeneracy between BR^a and 
the other parameters, in particular for the 10 GeV DM mass results for MIND, is mainly 
due to the 10 % systematic error assumed for the atmospheric background simulation and 
the energy resolution and statistics not being sufficient to tell the spectral shape of the 
background from a combination of BR„<7 and BR g <r. Notice that this effect also worsens 
the sensitivity estimated for the MIND detector, since relatively large values of BR„<r and 
BRq<7 simultaneously cannot be disentangled from a situation with a smaller atmospheric 
background accommodated within the 10 % systematic allowed. Thus, the sensitivity limits 
derived here, in particular for the MIND detector, are more conservative than the ones 
presented for Super-K in Ref. jlj] since only one branching ratio a time was considered there 
and in order for the signal to mimic the background a combination of both BR„<7 and a 
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FIG. IV: Typical results for the most probable regions of BR x cr (x = r, u, q) at 68 %, 90 %, and 
95 % posterior probability and a DM mass of 10 (upper pannels) and 25 GeV (lower panels). The 
solid lines show the GLACIER (thick) and 34 kt LArTPC (thin) results, while the dotted lines 
show the MIND results. The results correspond to 10 years of exposure for GLACIER, the 34 kt 
LArTPC and MIND. 



softer channel is necessary, as depicted by the degeneracy in the top left and top center 
panels of Fig. [TV] 



V. SUMMARY 



In this paper we have studied the detection of neutrinos as a possible indirect signal 
of low mass (10-25 GeV) DM annihilating in the Sun at the next generation of neutrino 
detectors, as well as the determination of which annihilation channels are present. Due to the 
results of the DAMA/LIBRA, CoGeNT and CRESST collaborations, this mass range, which 
is unaccessible to neutrino telescopes due to the low neutrino energy, is of great interest. 
We have focused mainly on the LArTPC and MIND detector technologies. Apart from 
having good energy and angular resolution capabilities, these two detector types are under 
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consideration for several future neutrino experiments. For the LArTPC, we consider both a 
34 kt and a 100 kt version, corresponding to detectors proposed for the LBNE and LAGUNA 
projects, respectively, and for the MIND we will consider a 100 kt detector, corresponding 
to the neutrino factory proposal or an upgraded version of INO. 



For the signal rate, we computed the so 
using the standard BP2000 Solar Model 



ar capture rate of DM following Ref. [l[ S] and 



77| . Our results are presented in the variables 



BR x cr (x = r, u, q), where BR X is the annihilation branching ratio into xx and a is the total 
spin-dependent DM-nucleon cross section. We have not included a systematic error on the 
signal, since this will mainly be an overall factor due to different capture mechanics and 
therefore lead to a simple rescaling of the BR x cr parameters. 

The main background to our signal is atmospheric neutrinos, for which we assume a 
systematic overall normalization error of 10 %. In order to suppress this background, we 
use the fact that the detectors under study have the capability of measuring the energy and 
direction of the incoming neutrino and demand that this direction is within the detector 
angular resolution from the direction of the Sun. 

Our results show that with 10 years data the facilities under study would have a sensitivity 
to BR X (T of (9(10~ 2 — 1) fb (c.f. Tab. [TTJ) with the best sensitivity being achieved for the 100 kt 
LArTPC detector and the worst for the MIND technology. This should be put into context 
by the current bounds from direct detection experiments, which give a < 400 fb at 10 GeV 
mass (PICASSO [15|) and a < 90 fb at 25 GeV (COUPP |lfj) and the region favoured 
by the DAMA /LIBRA signal around 10 GeV mass and a ~ 1 pb. The latest results 
by SIMPLE 17J improve on these limits down to ~ 10 fb. The sensitivities to different 
annihilation channels are essentially independent without correlations among them. For 
comparison with current experiments, the MIND detector would perform essentially on 
the level of the current Super-K js| data, although the analysis presented here is more 
conservative, allowing for degeneracies among the different channels and the atmospheric 
background that were not taken into account by the analysis of Super-K data and for which 
they would be potentially more dangerous, since only a counting experiment with no energy 
information was assumed in the study. On the other hand, the 100 kt LArTPC would 
out-perform Super-K by almost two orders of magnitude. In general, the sensitivities are 
better for the vv channel and worse for the qq channel, due to the spectra being harder and 
softer, respectively, as well as a lower event rate per annihilation for the qq channel, and 
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thus leading to different levels of confusion with the soft atmospheric spectrum. 

We have also simulated the results for a typical signal when BR x <t = 0.3 fb for all x. 
As expected from Tab. [TTJ it would be challenging for all of the considered experiments to 
distinguish BR g cx from zero, while the LArTPC detectors would clearly measure BR v u and 
to some extent BR r cr (c.f. Fig. lIVp . We point out that the experiments considered have 
the remarkable ability to tell apart the different annihilation branching ratios exploiting 
their energy resolution, something impossible with a counting experiment such as the study 
performed for present Super-K data. This demonstrates the complementarity of DM searches 
at neutrino detectors compared to other DM probes in the study of the nature of DM. 

We conclude that the prospects for improving the limits on low mass DM annihilations 
in the Sun in future neutrino experiments are in general good. Although current results for 
Super-K rule out the DAMA and CoGeNT signals with large branching ratios into channels 
giving rise to neutrinos, solutions with subdominant annihilation into these channels are still 
viable and present bounds can be improved by almost two orders of magnitude at a detector 
with the characteristics of GLACIER. Furthermore, these future experiments would have 
some capabilities of measuring the actual branching ratio into different channels by using the 
shape of the energy spectrum and such a measurement would provide important information 
on the nature of the DM. 
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